Hawking radiation and the boomerang behaviour of massive modes near a horizon 
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We discuss the behaviour of massive modes near a horizon based on a study of the dispersion 
relation and wave packet simulations of the Klein-Gordon equation. We point out an apparent 
paradox between two (in principle equivalent) pictures of black hole evaporation through Hawking 
radiation. In the picture in which the evaporation is due to the emission of positive-energy modes, 
one immediately obtains a threshold for the emission of massive particles. In the picture in which 
the evaporation is due to the absorption of negative-energy modes, such a threshold apparently does 
not exist. We resolve this paradox by tracing the evolution of the positive-energy massive modes 
with an energy below the threshold. These are seen to be emitted and move away from the black 
hole horizon, but they bounce back at a "red horizon" and are re-absorbed by the black hole, thus 
compensating exactly for the difference between the two pictures. For astrophysical black holes, 
the consequences are curious but do not affect the terrestrial constraints on observing Hawking 
radiation. For analogue gravity systems with massive modes, however, the consequences are crucial 
and rather surprising. 

PACS numbers: 04.70.-s, 04.70. Dy, 04.62. +v 



I. INTRODUCTION 

A. Black hole evaporation through the emission of 
positive-energy Hawking modes 

The evaporation of black holes through Hawking radia- 
tion is one of the cornerstones of post-classical gravity [I]- 
|3]. In Hawking's original derivation fl] for the case of a 
gravitational collapse, the emphasis lay on the positive- 
energy modes that cross the horizon just before it is ac- 
tually formed, escape from the black-hole spacetime and 
can in principle be detected by an asymptotic observer. 
Hawking already observes that the decrease of the black 
hole mass, and the accompanying decrease of the area 
of the event horizon "must, presumably, be caused by 
a flux of negative energy across the event horizon which 
balances the positive energy flux emitted to infinity. One 
might picture this negative energy flux in the following 
way. Just outside the event horizon there will be vir- 
tual pairs of particles, one with negative energy and one 
with positive energy. The negative particle . . . can tunnel 
through the event horizon to the region inside the black 
hole ... In this region the particle can exist as a real par- 
ticle with a timelike momentum vector even though its 
energy relative to infinity as measured by the time trans- 
lation Killing vector is negative. The other particle of 
the pair, having a positive energy, can escape to infinity 
where it constitutes a part of the thermal emission de- 
scribed above." Hawking warns, however, that it "should 
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be emphasized that these pictures of the mechanism re- 
sponsible for the thermal emission and area decrease are 
heuristic only and should not be taken too literally." 

With regard to massive particles. Hawking notes that 
"As [the black holes] got smaller, they would get hotter 
and so would radiate faster. As the temperature rose, 
it would exceed the rest mass of particles such as the 
electron and the muon and the black hole would begin to 
emit them also." Therefore, "the rate of particle emission 
in the asymptotic future . . . will again be that of a body 
with temperature Kjlix. The only difference from the 
zero rest mass case is that the frequency a; in the ther- 
mal factor (exp(27rajK~^) =F 1)~^ now includes the rest 
mass energy of the particle. Thus there will not be much 
emission of particles of rest mass m unless the tempera- 
ture njl-K is greater than m." A similar conclusion was 
reached in [3]. Indeed, for the black hole to emit for ex- 
ample an electron^ , its temperature must be on the order 
of T = 10^ Kelvin.^ Due to emission of massless particles 
the black hole will eventually become small enough for 
the temperature to reach 10^ K, and then the radiation 
will contain electrons and positrons. But for most of the 
lifetime of the black hole, the mass cutoff prevents any 
(significant) radiation of electrons. 

Implicit in the above reasoning is that the black hole 
emission corresponds (by definition) to what can be de- 
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tected by an asymptotic observer. This definition makes 
sense from a relativistic point of view because a black liole 
in relativity is, strictly speaking, not just the central ob- 
ject enclosed by a horizon, but also the entire spacetime 
surrounding it. This definition however has various well- 
known complications, for example that, strictly speak- 
ing, one needs an infinite amount of time to ascertain 
the existence of a black hole. In an astrophysical sense, 
therefore, the alternative definition in which the black 
hole is only what lies beyond the horizon makes much 
more sense. In the context of analogue gravity, which 
we will discuss in section |VII[ only the horizons are the 
essential ingredient and it makes little sense to speak of 
asymptotic observers. One can already anticipate that, 
from this second point of view, what is emitted by the 
black hole (i.e., what moves away from the horizon) is 
not necessarily equal to what is detected at infinity. As 
we will see, the difference can actually be substantial. 



B. Black hole evaporation through the absorption 
of negative-energy Hawking modes 

Hawking's heuristic picture in which the black hole 
shrinks through the absorption of negative-energy modes 
was given a physical basis in the following context. 
An important problem with Hawking's calculation is 
the so-called "transplanckian problem" : the exponential 
frequency-shift near the horizon means that particles de- 
tected just seconds after the formation of the black hole 
seemingly originate from modes with an energy exceed- 
ing by far the total energy of the observable universe 
just before the formation of the black hole. In 1981, in 
what is usually considered the founding paper of analogue 
gravity (see for a review), Unruh |5] noted that the 
propagation of sound waves in a fluid with a subsonic- 
to-supersonic transition is formally identical to the prop- 
agation of scalar fields in a black hole spacetime. Since 
the continuum description of fluid systems has a natural 
cut-off (the intermolecular distance), the transplanckian 
problem of Hawking radiation can be studied in this fluid- 
analogy context. Unruh took up his own suggestion and 
studied the (typically "subluminal" ) modification of the 
dispersion relation for sound waves compared to a "rela- 
tivistic" or acoustic one, and its influence on the Hawking 
spectrum f^. For such subluminal dispersion, an outgo- 
ing positive-energy Hawking packet can be traced back to 
a mixed positive/negative-energy mode originating from 
outside the horizon. If the dispersion is changed to su- 
perluminal then the positive-energy mode originates from 
inside the horizon, but this difference does not alter the 
size of the Hawking effect compared to the subluminal 
case. This suggests that the evaporation (loss of mass) 
of the black hole does not rely so much on this positive- 
energy emission, but instead has its physical origin in the 
absorption of the negative-energy partner. The negative- 
energy mode is always consumed by the black hole, re- 
gardless of the details of the high-energy dispersion. 



It is far beyond our intention to clear up the exact 
physical origin of Hawking emission and black-hole evap- 
oration. What we do want to point out, however, is 
the following. The emission of massive positive-energy 
modes immediately seems to be limited by a lower cut- 
off, as pointed out above, since the energy of a massive 
mode in the flat space-time of an asymptotic observer 
cannot be less than its mass-energy. This is not the 
case, however, for the absorption of massive negative- 
energy modes, since these fall into the black hole and 
so never obey the flat space-time dispersion relation. In 
fact, as we will discuss in detail, the process of creation 
of positive/negative-energy pairs near the horizon is not 
restricted by the mass of the corresponding modes. Mas- 
sive negative-energy modes are created and absorbed by 
the black hole according to exactly the same mechanism 
and at exactly the same rate'^ as massless modes. Since 
the two pictures of black hole evaporation through Hawk- 
ing radiation (either through the emission of positive- 
energy modes or through the absorption of negative- 
energy modes) should in principle lead to exactly the 
same evaporation rate, this raises an apparent paradox. 
The paradox is qualitatively easy to solve: all positive- 
energy modes emitted with energy (relative to asymp- 
totic infinity) below the rest-mass energy must somehow 
be reabsorbed by the black hole in order to counterbal- 
ance the "excess" negative-energy absorption. A detailed 
analysis of how this process takes place leads to some 
curious consequences, and in particular to several poten- 
tially crucial aspects for experiments in analogue gravity 
systems with massive modes. 

Before starting our analysis, note that we will talk 
about massive modes to avoid any caveats concerning 
the definition of particles in a curved spacetime. This is 
not just a matter of terminology. Strictly speaking, par- 
ticles do not go on-shell until they have reached the flat 
region of spacetime (either at infinity, or when they are 
trapped by the potential well of a particle detector before 
reaching infinity). Care should therefore be taken when 
reasoning in terms of particles. For example, the usual 
heuristic argument about particle pair creation near the 
horizon (due to the effect of curvature on a distance of the 
order of the Compton wavelength) might induce one to 
think that massive particle pairs are created near a black 
hole horizon only when the action is sufficient to provide 
an energy at least equal to the pair's mass. This argu- 
ment might be useful in that it correctly reproduces the 
threshold for particles to appear at infinity, as mentioned 
above. However, in order to make estimates about what 
happens anywhere between the black hole horizon and in- 
finity, it can be quite deceptive, and it is more indicated 
to reason in terms of modes or wave packets, as we will 
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do. For analogue gravity systems, such complications do 
not usually arise because the emitted modes are in gen- 
eral phonons or other collective excitations, which are 
in any case naturally decribed in terms of waves rather 
than particles. Moreover, the configurations considered 
for experimental detection of analogue Hawking radiation 
are typically stationary and have an asymptotically con- 
stant background flow velocity U. The ambiguity of the 
quantum vacuum definition in a non-stationary curved 
spacetime [2] does therefore not arise here (see [6] for 
useful remarks in the case of phonons in a BEC) . We will 
come back extensively to the question of detectability in 
sections IVll and IVni 



II. DISPERSION RELATION FOR MASSIVE 
MODES IN A BLACK HOLE SPACE-TIME 

The general dispersion relation for massive Klein- 
Gordon modes trying to escape a black hole can be writ- 
ten 



lOJ- 



Ukf 



m 



,2;„2 



c'k 



(1) 



Note that, strictly speaking, the mass M of the mode 
is given by M = mh/c^. We set h = 1 and will refer 
throughout to m as the mass. As usual, U < can 
either represent the velocity of an observer freely falling 
into the black hole (in the case of gravity) , or the velocity 
of a background flow against which the mode counter- 
propagates (for analogue gravity systems). Both cases 
are described by exactly the same Painleve-Gullstrand- 
Lemaitre metric 



[c^ - U{x)^]dt^ - 2U{x)dtdx - dx^ 



(2) 



where x is the direction perpendicular to the horizon 
(which contains all the essential features that we are 
interested in**), and c represents the speed of light, or 
in analogue gravity the (low-frequency) speed of sound 
or any other characteristic velocity of the system which 
leads to an effective relativistic behaviour (typically) at 
low frequencies. The strict equivalence (at least at low 
frequencies) means that we do not have to worry about 
terminology, and so we will typically refer quite gener- 
ally to, for instance, \U\ = c as the "sonic" case, U as the 
"background flow" , and "counter-propagating modes" for 
modes that try to escape the black hole (or enter a white 
hole), as if they were moving against a background flow. 



* Note that a bouncing effect similar to what we will describe also 
exists for (even massless) modes with non-zero angular momen- 
tum, see e.g. |10| . as is well known. The crucial difference is that 
the bouncing which we will describe is valid even for modes that 
move away from the black hole in a purely radial fashion. It is 
therefore truly a consequence of the mass of the modes. 



The action of the Klein- Gordon field 
(U) is 



in the metric 



S = 



dt dx 



which gives the field equation 



(3) 



inn?(j) = 
(4) 

and the dispersion relation ([!]) . Wave packets obeying the 
Klein-Gordon equation Q possess two conserved quan- 
tities due to the invariance of the action ([3| under 1) the 
transformation (f) — >■ e*"(/), a constant, and 2) time trans- 
lation (for time-independent U). The former invariance 
gives conservation of the Klein-Gordon norm 



N = 



2c2 



(5) 

whereas the latter gives conservation of (pseudo-)energy 



E = - 
2 



dx 



(6) 

For waves packets confined to a region where the flow 
velocity U is constant, the Klein-Gordon norm ([S]) can 
be written in fc-space in terms of the Fourier transform 
4){k) as 



N : 



1 r°° 

-/ dk{ij-uk)\m\^ (7) 



while in similar circumstances the energy (|6| takes the 
form 



E 



dkcoiuj -Uk)\4){k)(' 



(8) 



The quantity that determines the sign of the norm ([s]) 
of a wave is the sign of its frequency uj — Uk in a frame 
co-moving with the background flow U. For waves with 
positive values of the frequency w in the black-hole/lab 
frame, the sign of the (pseudo-) energy ^ is also given 
by the sign of the co-moving frequency lo — Uk. 

A key point is that the dispersion relation M far away 
from the black hole (in Minkowski space-time) is 



UJ 



m + ck 



.2;.2 



(9) 



which does not allow any particles with cj < m. However, 
sufficiently close to the black hole horizon U — >■ — c, so 
from ([ij massive modes can in principle be created even 
at frequencies < w < m. 

We will now plot the dispersion relation ([I]), i.e. 



•-Uk = ±y^, 



(10) 



for values of a; > greater than or less than the mass 
m. For fixed uj the left- and right-hand sides of ( 10 ) are 
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(a) 



Lo-Uk 





uj > m 



uj < m 





= c 



u) > m 



UJ < m 



FIG. 1: Positive-frequency part of the dispersion relation ( 10 1 
for u) > m and (a) 17 = 0, (b) various values of 17 < 0. 
fc < corresponds to co-propagating modes, fc > to counter- 
propagating ones. The massless sonic branches uj — Uk = ±cfc 
are indicated (dashed red lines) for comparison. Qualitatively, 
the u) > m massive case is identical to the massless case. 



plotted as functions of fc; the intersection points of these 
plots are the solutions of the dispersion relation. Fig- 
ures [l] and [2] represent the positive-frequency branch of 



the dispersion relation (plus sign in ( 10 )), while Fig. 
picts the negative-frequency branch (minus sign in (10|). 
Note that "positive/negative- frequency" always refers to 
the comoving frequency (uj — Uk), and therefore coincides 
with the concept of positive/negative norm, as just men- 
tioned. We always draw U < {a flow moving to the left) 
and will focus on the modes counter propagating against 
the flow, i.e. modes moving to the right in a frame co- 
moving with the flow. The counter-propagating modes 
have positive phase and group velocities in the frame co- 
moving with the flow: {uj — Uk)/k > 0, d{uj — Uk)/dk > 0. 
Our interest is in counter-propagating modes because 
these modes try to escape the black hole (or enter a white 
hole). 



FIG. 2: Positive-frequency part of the dispersion relation 1 10 1 
for < UJ < m and various values oi U < 0. For the value 
of |f/sub| < c shown in (a), there are two counter-propagating 
(fe > 0) modes, one with a positive group velocity duj/dk > 
in the black-hole/lab frame, the other with a negative group 
velocity in this frame. In (b), U^nh = U* , the critical value at 
which these two solutions merge and only one (double) root k* 
remains. For |[/sub| < \U*\, there are no longer any positive- 
frequency modes. U* is therefore characteristic of a turning 
point or "red horizon" (see main text). The massless sonic 
branches u — Uk — ±cfe are again indicated (dashed red lines) 
for comparison. Note that all co-propagating modes {k < 0) 
have disappeared, while the number of counter-propagating 
modes (k > 0) depends on the value of U. 



A. ij > m (positive-frequency part) 

The case w > m is very similar to the standard mass- 
less case, see Fig. [T] In the absence of a counter-flow, i.e. 
far away from the black hole (Fig. [ija)), there are two 
positive-frequency roots ki and k2 = —ki with the same 
modulus and opposite direction of propagation. When 
U =^ (Fig. fljb)), the situation depends on the value of 
U: for |C/subl< c, there are two positive- frequency roots 



fci < and k2 ^ -ki, for |J7sonic| = c and |J7supcr| > c, 
only one positive-frequency root remains, namely the 
mode fci < co-propagating with the background flow 
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w - Uk 

I Usonic I C ' < 




FIG. 3: Negative-frequency part of the dispersion rela- 



tion ( 10 1 for oj > 0. The behaviour is independent of whether 



Lj < m or Lo > m, and is qualitatively identical to the mass- 
less case (in dashed red line). For IL^subl < c and IL^sonicI = c, 
there are no negative-frequency solutions. A single counter- 
propagating negative-frequency root appears for | (/super] > c. 



U. For I [/super I > c, the second root k2 has actually been 
converted into a counter-propagating negative-frequency 
root, see Fig.|3] This conversion from positive to negative 
frequency across the horizon is precisely the essence of 
the Hawking mechanism (at least in the case where there 
are no higher-order fc-terms in the dispersion relation). 
For |C/sonic| = c, ^2 has "disappeared": /c2 — > +oo when 
\U\ — > c^, and reappears as fc2 ~^ ~oo when \U\ — > c+. 
This infinite frequency-shift at the horizon (|C/sonic| = c) 
is characteristic of the transplanckian problem, which 
(just like in the massless relativistic case) is present in our 
problem since we do not consider higher-order fc-terms 
in the dispersion relation in our present discussion. Note 
that the exact values of the different roots will be slightly 
different compared to the massless case, but in qualita- 
tive terms the behaviour for the positive-frequency modes 
associated with a; > m is identical to the massless case. 



ground flow and has negative group velocity relative to 
the black hole. For a critical value U* of U (Fig. [2]^b)), 
these two solutions merge and only one (double) root 
k* remains. For even smaller values of \U\, as in the 
case oi U ~ 0, there are no longer any positive-frequency 
solutions. These two fc > roots that merge at k* cor- 
respond to the same counter-propagating mode; it prop- 
agates away from the horizon (positive group velocity), 
experiences a turning point when the flow has decreased 
to U* , and is dragged back into the black hole (nega- 
tive group velocity). The positive- frequency mode with 
< a; < TO thus exhibits a boomerang trajectory: after 
being emitted away from the horizon it returns, and falls 
into the black hole. This boomerang behaviour is fur- 
ther illustrated below by a ray plot and a wave packet 
simulation. 



C. Negative- frequency part (all values of to) 

The negative-frequency part is shown in Fig. [3] The 
behaviour of the negative-frequency modes is indepen- 
dent of the value of w > 0. It is essentially the 
same as in the massless case: inside the horizon (where 
\U\ = I [/super I > c), a single negative- frequency solu- 
tion exists. Outside the horizon (|[/| < c), there are 
no negative-frequency solutions. This corresponds pre- 
cisely to the "disappearance" of the positive-frequency 
fc2-solution in the massless and lo > m cases: k2 converts 
from a positive-frequency root into a negative-frequency 
root across the horizon. 

The key point, however, is that this same appearance 
of a negative-frequency root when \U\ > c is also valid for 
the case co < m. Moreover, since this negative-frequency 
mode is absorbed by the black hole and hence moves into 
a region where \U\ increases, its existence is allowed at 
all times during its subsequent evolution. 



III. DETERMINATION OF U* 



B. u) < m (positive-frequency part) 

For Lu < m, important differences with the massless 
case show up, see Fig. [2] For U ~ (not shown), there 
are no positive-frequency uj < m solutions to the disper- 
sion relation, in accordance with For |[/sonie| = c and 
|C4upor| > c, there is only one positive- frequency root. 
Curiously, and contrarily to the case w > m, the single re- 
maining mode here is a counter-propagating mode, while 
the mode co-propagating with the background flow is for- 
bidden. For I [/sub I < c, the number of roots depends on 
the value of U: for sufficiently large \U\ (Fig.[2]ja)), there 
are two positive-frequency roots. These are both counter- 
propagating, but only one of them (with the higher k) has 
a positive group velocity du>/dk > in the black-hole/lab 
frame, while the other one is dragged along by the back- 



We have seen graphically that in the case w < to, 
there exists a critical background flow velocity U* such 
that, for \U\ > \U*\, there are two (counter-propagating) 
positive-frequency solutions, whereas for |[/| < \U*\, 
both these solutions disappear. For \U\ — \U*\, there 
is therefore a double root k* or a saddle-node bifurcation 
in dynamical-systems language [TTl [T2] . 

For 7^ , the dispersion relation ^ can be written 

as 

, 2ujU lo"^ ~ m? , , 

^ -[7^^+T7^-«' (11) 

with a discriminant 

^-([7^3^-4([7^^(^ (^2) 
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There will be real solutions for k as long as A > 0, i.e. 



C/2 >c2 



1 — {uj/mY 



ki 2 is then 



fcl 2 = 



t/2 



The general solution for the roots 



(a;2-m2)([/2-c2) 
w2i72 



1± \ 1 



(13) 

The critical value U* can be obtained by requiring A to 
vanish. This can only be satisfied when a; < m, and gives 



U*=±csll-[-) , (14) 



with a corresponding critical wavenumber 



1 - 



(15) 



As we consider flows U{x) < 0, the lower signs apply in 
(14) and (15). The critical value (14) occurs outside the 



black-hole horizon U = —c and corresponds to a turning 
point for the counter-propagating lo < m mode as it tries, 
and fails, to escape the pull of the black hole. 

The case U = ±c gives in the dispersion relation ([I]) 



m 



2ujc 



(16) 



so that only one solution remains (for U < 0). The sign 
of the single root depends on the relation between uj and 
TO. Taking U — J/sonic = — c, the solution (16) has k < 
for w > TO, corresponding to a co-propagating mode 
crossing the horizon into the black hole. For w < to. 



k > in (16), corresponding to a counter-propagating 



mode being dragged across the horizon after encountering 
a turning point while trying to escape to infinity. These 
quantitative results are in agreement with Figs, [ijb) and 



m 



If w = 



then 



2mU 



[/2 



(17) 



and the case is qualitatively the same as a; > to, except 
that there is always one solution at fc = 0. 



FIG. 4: Ray trajectories for < to < m (top) and u > 
m (bottom) modes trying to escape from a black hole on 
both sides of the horizon (green: outside the horizon, blue: 
inside the horizon). After remaining extremely close to the 
horizon for some time, the modes that started off just inside 
the horizon (blue lines) are dragged into the interior of the 
black hole. The modes just outside the horizon eventually 
manage to escape at least some distance from the black hole. 
For the case uj > m (green ray on the bottom), they escape 
to infinity. For uj < m (green ray on the top), the mode 
escapes some distance from the horizon, but reaches a turning 
point or "red horizon" , bounces back and is re-absorbed. The 
ray trajectories can also be interpreted as corresponding to 
modes trying to penetrate a white hole, by inversing the time 
direction. The fiow then reverses direction {U > 0) and the 
interior of the white hole is on the left of the figures. The 
counter-propagating rays pile up on the horizon from both 
sides (green: from the outside, blue: from the inside). Note 
that, for the case ui < m (top), the green ray originates inside 
the white hole, is ejected, turns around at the "blue horizon" 
(the time-reverse of the red horizon, see also main text), tries 
to enter the white hole again but is unable to penetrate across 
the white-hole horizon where it piles up. 



IV. RAY PLOTS AND WAVE PACKET 
SIMULATION 

We will now plot ray trajectories for the cases u < m 
and UJ > m. This will confirm that, for ut < m, there is 
a (mass- and frequency-dependent) turning point, where 
the outgoing positive-frequency mode is blocked, bounces 
back towards the horizon and is reabsorbed by the black 
hole. 



A. Ray trajectories 

We plot in Fig. |4] ray trajectories for counter- 
propagating modes that start off infinitesimally close to 
the horizon on both sides of it. The mode starting off on 
the inside eventually falls down into the interior of the 
black hole, independently of its frequency w, since |C/| > c 
inside the horizon. For modes starting just outside the 
horizon, the behaviour is completely different depending 
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on whether uj > m ov uj < m. 

For a; > m, as in the massless case, modes starting 
off just outside the horizon eventually escape to infin- 
ity. Note that the group velocity of the massive mode 
is slightly smaller than that of its massless counterpart, 
and the difference increases towards low k, as can be seen 
graphically from the slope of w — Uk in Fig. [T] How- 
ever, on the outside of the horizon, \U\ < c, and so 
duj/dk > \U\, as can again easily be confirmed graphi- 
cally from Fig. [T] This can be interpreted as follows: the 
mode starts off with a very high wavenumber, at which 
the group velocity dio' /dk in a frame comoving with the 
flow [uj' = oj — Uk) barely differs from the massless case, 
i.e. duj'/dk{x ~ xh + e) ~ c > U{x ^ xh + e), with 
xh the position of the horizon and e a small positive 
quantity. The mode therefore gradually escapes, and is 
redshifted (towards lower k) in the process. The decrease 
in \U\ as the ray moves away from the horizon is more 
important than the decrease in duj' /dk (graphically: the 
slope duj' /dk is always larger than \U\), and so the mas- 
sive mode indeed escapes. 

For w < m, on the other hand, as the massive mode 
ki moves away from the horizon and is redshifted in the 
process, it will reach the critical U* beyond which there 
are no longer any positive-energy solutions for k. The 
process here is typical for a saddle-node bifurcation: at 
U* the two roots ki (with positive group velocity Cg — 
duj/dk) and k2 (with negative Cg) merge and Cg vanishes. 
The outgoing ray is therefore blocked and turns back 
towards the black hole, in a boomerang fashion. In other 
words, the outgoing (cg > 0) mode ki is converted (red- 
shifted) into an ingoing (cg < 0) mode ^2 when it reaches 
the point where [/ = [/*, and is subsequently re-absorbed 
by the black hole. Note that the counter-propagating 
character of the mode with respect to a background flow 
is essential for this behaviour, just like a boomerang must 
be thrown against the wind direction (typically at 45°) 
in order to come back. 

We call the turning point where U = U* a "red hori- 
zon" in line with the term "blue horizon" introduced 
in |12j : the red horizon is a turning point for red-shifted 
waves, and in the whole process at the red horizon the 
wave undergoes further redshifting. We could also re- 
verse the time direction in Fig. |4]to study what happens 
for a mode trying to enter a white hole. In that case, the 
blocking line or turning point would be associated with a 
blue-shifting process, and we would recover a blue hori- 
zon in precisely the same sense as introduced in [12 . A 
red horizon in the black-hole case thus corresponds (by 
time-inversion) to a blue horizon in the white-hole case. 
Note that the location of such a red/blue horizon is auto- 
matically frequency-dependent, as we discuss below, even 
though no higher-order terms in k have been introduced 
in the dispersion relation. 

The bouncing behaviour of the green ray in the top 
plot in Fig. [4] (0 < w < m) , together with the continual 
red-shifting in time, is shown by a numerical solution 
for a Klein-Gordon wave packet centered on this ray, see 



Fig. [Sj The wave packet at various points in time is 
shown in Fig. [5j with time increasing from top to bottom. 
Initially the wave packet is outside the horizon (the grey 
vertical line) and moves to the right, away from the black 
hole. The packet reaches a maximum distance from the 
horizon before reversing direction and falling through the 
horizon into the black hole. The extreme red-shifting of 
the wave packet that accompanies its motion is evident. 
(The numerical simulation was performed backwards in 
time, beginning with the packet inside the black hole and 
solving for its history; equivalently, the simulation was 
performed forwards in time for the white hole obtained 
by reversing the flow U .) 



B. Change of red horizon with changing mass or 
frequency 

The evolution of the red horizon with mass m and fre- 
quency u! is shown in Fig. |6] Note that the existence of 
the red horizon assumes < w < m. 

When the mass m is increased, the position of the red 
horizon moves closer to the black-hole horizon. Con- 
versely, if m is decreased (keeping m > uj), the position of 
the red horizon moves away from the black-hole horizon. 
This is easy to see from the value of the critical U* which 
identifles the red horizon, see eq. (14): as m increases, 
\U*\ increases, and so lies closer to^I/| = c, i.e. to the 
black-hole horizon. A massive mode with m — > 00 will 
not be able to escape any distance from the horizon, since 
|C/*| — > c. In the lower limit to — > w, on the other hand, 
we obtain U* — > 0, i.e. the mode can asymptotically 
reach infinity. 

If UJ is increased (keeping uj < to) , the position of the 
red horizon moves further from the black-hole horizon. 
Conversely, if uj is decreased, the position of the red hori- 
zon moves closer to the black-hole horizon. This is again 
easy to see from eq. (14). In the limit w — >■ 0, the mode 
cannot escape any distance from the horizon (|C/*| — )■ c), 
whereas in the limit w — to, we obtain |C/*| — >■ 0, so that 



the mode escapes to infinity. 



V. PHASE DIAGRAM 

Some of the preceding considerations can be sum- 
marised in a phase diagram of fiow speed \U\ and pe- 
riod T = 2n /uj, see Fig. [Tj The purple horizontal dashed 
line in Fig. [t] is \U\ = c, which represents the black- 
hole horizon. It also represents a "negative horizon" |T2] : 
negative-frequency modes can exist only for \ U\ > c. This 
negative horizon is indicated by a horizontal black line in 
Fig. [7j although the negative horizon coincides with the 
black-hole horizon in the case discussed here, these two 
horizons can separate, for example in the case of surface 
waves on moving water |12) . 

For massless modes, the purple (and black) line 
would strictly mark the only blocking line for modes: 
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FIG. 5: Evolution of a Klein-Gordon wave packet with frequency spectrum contained in the range < w < m. The grey 
vertical line is the horizon of a black hole (interior on the left, exterior on the right). The packet starts (top) close to the 
horizon, propagates away from the black hole, bounces and falls into the black hole interior, red shifting all the way. 
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FIG. 6: Mass- and frequency-dependence of the "red hori- 
zon". Mass (top): mrod > m-grocn > TTiyoiiow The heavier the 
mode, the less distance it can escape the black hole before 
being bounced back and reabsorbed. Frequency (bottom); 
ojrcd > i^grcon > t^yoiiow . The higher the frequency, the further 
towards asymptotic Minkowski space the mode can escape. 
The green curve in both graphics corresponds to the values 
used in Fig. [4] 

all counter-propagating massless modes have a positive 
group velocity duj/dk > in the lab frame as long as 
\U\ < c (they overcome the background flow), but a 
negative duj/dk for |C/| > c (they are dragged along by 
the background flovif). |J7| = c therefore marks the only 
blocking line {duj/dk = 0) for modes in the massless case. 

For massive modes, the picture becomes more compli- 
cated. The red line in Fig. [7] shows the critical speed 
\U\ = \U*\ at the red horizon as a function of T. The 
vertical dashed green line corresponds to w = m, i.e. 
T = Tm = 2iT/m. As T approaches T,„ from above the 
red horizon moves further away from the black-hole hori- 
zon |[/| = c and reaches spatial infinity at T = Tm, in 
line with the results of the last section. For T < r„j, 
the red horizon has disappeared and the behaviour is 
qualitatively the same as the massless case: counter- 
propagating modes have a positive lab-frame group ve- 
locity doj/dk > when \U\ < c, but are dragged along 
by the background flow {duj/dk < 0) when \U\ > c. The 
only blocking line {dui/dk — 0) for T < Tm is there- 
fore \U\ = c. For T > Tm, however, a second blocking 
line \U\ = \U*\ is present, the "red horizon": for a given 
period T (frequency w), massive modes can only exist 
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FIG. 7: Phase diagram of flow speed \U\ versus period 
T — 2tt/uu. The purple horizontal dashed line is \U\ — c and 
corresponds to the black hole horizon. The red line shows the 
critical speed \U*\, which defines the red horizon. The verti- 
cal dashed green line corresponds to a; = m, i.e. Tm = 27r/m. 
For T < Tm, everything still behaves qualitatively as in the 
massless case. For T > Tm, a second turning point of zero 
group velocity appears, the red horizon. Massive modes can 
only exist for \U\ > \U*\. The numbers in the phase diagram 
represent the number of roots (fc's) in each region and their 
signs; the upper position for a sign means the root has posi- 
tive (co-moving) frequency, the lower position means the root 
has negative frequency. 



for \U\ > \U*\, i.e. when the counterflow is sufficiently 
strong. For weaker background flows \U\ < \U*\, their 
existence is prohibited by the dispersion relation ([T]). As 
the mode moves away from the horizon (\U\ decreases 
below c), it is blocked {duj/dk = 0) at the red horizon 
\U\ = \U*\, bounces back {\U\ increases again), crosses 
the black- hole horizon {\U*\ = c) and is absorbed by the 
black hole. 

The numbers on the phase diagram represent the num- 
ber of roots (values of k) and their signs; the positions of 
the signs show whether the root is positive (co-moving) 
frequency (upper position) or negative frequency (lower 
position). For example, a mode with T > Tm has roots 
for < |C/*|, one (double) root k > Q iov U = U* and 
two positive roots (fci, k2 > 0) for |C/*| < |C/| < c. When 
|[/| — c~, fci — > -|-oo, and reappears as fci — > — cx) when 
\U\ — ?► . For \U\ > c, there is therefore one negative 
and one positive root. For T < Tm, there are 2 roots, 
one positive {ki, counter-propagating) and one negative 
(/c2j co-propagating), both with positive (co-moving) fre- 
quency, as long as \U\ < c. For \U\ > c, the fci solution 
has converted into a negative-frequency and negative-/c 
root and therefore ki,k2 < 0. As before, the transi- 
tion from /ci > to /ci < goes through ki ~> ±oo for 
\U\ — )■ c^, i.e.: fci effectively disappears for \U\ = c. 
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VI. HAWKING RADIATION 

As was seen in the foregoing, the behaviour of aU mas- 
sive (counter-propagating) modes, including those with 
w < m, is identical very close to the black-hole hori- 
zon: they all experience the familiar infinite blue-shifting. 
This behaviour at the horizon is in fact a serious prob- 
lem for the derivation of the Hawking effect for real 
black holes, as it invalidates the assumption of negligi- 
ble back-reaction on the black hole. Nevertheless it is 
now well-established that the introduction of non-linear 
dispersion at high wave numbers does not destroy the 
Hawking mechanism (although see [131 [H] for some pos- 
sible complications with respect to the vacuum selection 
in the case of super-luminal dispersion), while it limits 
the blue-shifting at the horizon. This removal of infinite 
(co-moving) frequencies through dispersion is certainly 
present in all analogue systems and gives a robust basis 
for investigations of Hawking radiation. As all massive 
modes are not distinguished by their behaviour at the 
horizon, where the Hawking mechanism occurs, any fields 
described by the Klein-Gordon equation will exhibit the 
Hawking effect at all frequencies u, i.e. without a cutoff 
given by the mass. 

But the fact that a particular mode exhibits the Hawk- 
ing effect does not imply that it contributes to the evap- 
oration of the black hole. After being emitted from the 
black hole, particles with ui < m eventually turn around 
and fall though the horizon, so they do not contribute 
to the evaporation. But one could imagine capturing the 
u! < m particles before they fall back into the black hole, 
in which case they would contribute to the shrinking of 
the horizon. A detector D placed sufficiently close to 
the black hole will detect a significant amount of massive 
particles emitted by the black hole that are not detected 
by a detector at infinity. The same detector D would 
also register massive Hawking particles falling back into 
the black hole. 

For a Schwarzschild black hole, the position x* of the 
red horizon is obtained from {U* /c)'^ = Rs/x*, with Rs 
the Schwarzschild radius and U* given by eq.(14|. This 
leads, for example, to x* « 3i?s for uj = 0.8m, and 
X* « 50i?s for UJ — 0.99to. The spectral resolution of 
a hypothetical particle detector placed at a position x^ 
will essentially be determined by the fiy-by time of the 
emitted mode: the time between crossing the detector on 
its way out and its disappearance behind the black hole 
horizon, after having bounced on the red horizon. To ob- 
tain a sufficient resolution, one would need a supermas- 
sive black hole, and/or place the detector extremely close 
to the black hole horizon, where the mode is slowed down 
by the black hole's gravitational pull. In that region, all 
sorts of additional astrophysical effects would probably 
put strong limits on the detectability of the boomerang 
behaviour. Although these issues are not very different 
from the general problems of detecting (massless) Hawk- 
ing radiation, it is clear that such considerations are of 
theoretical interest only in the astrophysical case. 



However, for an analogue system, there is no issue with 
the spectral resolution since the modes can perfectly well 
be detected inside the horizon. For an analogue system in 
which a field obeys the Klein-Gordon equation, the lack 
of a mass cutoff in the emission of Hawking radiation is 
very important, and we will dedicate the next section to 
it. 

Explicit quantitative estimates for the detection rate 
W of massive particles with ksTn ^ lu < m (with Th 
the usual Hawking temperature) can be obtained through 
the semiclassical tunneling formalism |151 116j . This gives 
the standard Hawking result for massless particles, mod- 
ulated by a factor which depends on the particle mass 
m, the difference m? — uj'^ , the distance to the detec- 
tor, and the black-hole spacetime profile U{x). While the 
general expression is rather involved, it strongly simpli- 
fies when xd ^ xh, with xh the radius of the horizon. 
In that case, one obtains the simple approximation 









oc exp 




exp 



2xd_ 
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(18) 



(in units such that /i = 1, as before), with Th ~ 
\dU / dx\x=XH 1'^''^ tti^ standard Hawking temperature. 
Compared to the standard case, the probability that a 
massive particle with uj < m arrives at a distance xd 
decreases exponentially with xd, as expected. 



VII. ANALOGUE GRAVITY 

The study of Hawking radiation in analogue-gravity 
systems opens up possibilities that are excluded in as- 
trophysics. One can consider black- and white-hole hori- 
zons, which both produce Hawking radiation, and indeed 
combinations of horizons that communicate with each 
other (for example, the so-called black- hole laser |17)). 
One can also perform measurements on both sides of a 
horizon, so detection of both particles in a Hawking pair 
and measurement of their correlation becomes possible. 

The study of astrophysical black holes may lead to the 
simplistic view that there is a lower frequency (energy) 
threshold in the production of massive particles by hori- 
zons. This in turn would lead to the view that an ana- 
logue system where a field described by the Klein- Gordon 
equation experiences a horizon is a poor choice for observ- 
ing the Hawking effect. In reality the Hawking radiation 
produced by the horizon is in no way suppressed by the 
mass, and such an analogue system is in this regard just 
as promising a candidate as one in which the field obeys 
the massless wave equation. The u < m particles would 
be just as accessible to detection as the uj > m particles. 
In the case of a black-hole horizon the uj < m particles 
end up going into the "black hole", whereas the uj > m 
particles propagate away from the other side of the hori- 
zon, but all can be detected. For a white-hole horizon, 
what happens depends on how the dispersion relation is 
modified at high k in the particular analogue system. All 
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real systems will be dispersive at the high k values expe- 
rienced by modes near the horizon and instead of sticking 
at the white-hole horizon the rays will drift away from the 
horizon on one side or the other depending on whether 
the dispersion is super-luminal or sub-luminal at high k. 
Again, there is no problem in principle in detecting all 
the particles, for all uj, in this white-hole case. 

One may wonder whether the boomerang effect will 
not be mixed with the — possibly very complicated — 
dispersive high-fc behaviour near the horizon. Two ele- 
ments should be borne in mind. The first is that detailed 
studies based on rather general (mainly "subluminal" , 
see e.g. [9]) dispersion relations as well as for concrete 
systems with massless waves (BECs, surface waves, etc.) 
have shown that such potential ultraviolet problems do 
not spoil the Hawking effect and a thermal spectrum can 
be recovered under quite general conditions on the veloc- 
ity profiles. For BECs, for example, this is demonstrated 
starting from the Bogoliubov-de Gennes equations ^6j, 
and the outcome can safely be expected to be indepen- 
dent of the masslessness or massiveness of the phonons. 
The second key point is that the boomerang effect is de- 
tectable as long as the background flow velocity U is non- 
zero. In analogue systems, this region can be extended 
to a distance far beyond the immediate surroundings of 
the horizon, so that any complicated high-A: behaviour 
has been completely smoothened out. 

It remains to identify analogue-gravity systems where 
these ideas could be tested experimentally. In most 
analogue systems that have been considered the waves 
obey the massless wave equation at low wave-numbers 
(with modifications at high wave- numbers) . An excep- 
tion is a recent study of acoustic waves in a rotating ion 
ring [IHl [H] where the discreteness of the ions leads to 
a dispersion relation at low k that is not of the mass- 
less form, but these waves are not Klein-Gordon waves. 
We give three examples of systems where waves obey the 
Klein-Gordon equation and where the creation of hori- 
zons for the waves is possible. 

Bose-Einstein condensates (BECs) are among the most 
discussed systems for analogue gravity. Theoretically, 
they provide the advantage of being well-studied and 
rather well-understood, while experimentally they are 
clean and quite flexible. Black hole configurations in a 
BEG were recently realized for the first time |2Q]. The 
phonons that are obtained in standard BEGs are mass- 
less, but in 21J a scheme was described for producing a 
"massive" phonon in a two-component BEG. In the pres- 
ence of a horizon created by a sub- to super-sonic fiow 
these massive phonons should then behave as described 
above (with an appropriate modification at high fc). 

The second example is Langmuir waves in a moving 
plasma [55]. These waves have a dispersion relation of 
the form (w — Uk)"^ = uip +ujpRjjk'^, with ujp the plasma 
frequency and Rd the Debye radius. This is the Klein- 
Gordon dispersion relation with the cut-off frequency Up 
playing the role of mass. One obvious undesirable feature 
in this system is Landau damping [221 of the waves, but 



all analogue-gravity systems present their own challenges. 

The third example is from barotropic waves [55] in 
fiuid mechanics. Barotropic waves are an assortment of 
waves in fluids with rotation, one type of which, known as 
inertia-gravity or Poincare waves, has the Klein-Gordon 
dispersion relation [23 = + ghk^, where / is the 
Goriolis parameter, g the gravitational acceleration and h 
is the fluid depth. As in the previous examples, a horizon 
can be created by a fluid flow. 

The presence of a lower frequency cutoff in the flat (but 
not curved) space-time dispersion relation is the feature 
of Klein-Gordon waves that leads to the interesting be- 
haviour we have discussed in this paper. This lower cutoff 
is not unique to Klein- Gordon waves and we note an- 
other possible analogue-gravity system with waves that 
exhibit a lower frequency cutoff, namely one based on 
spin waves. Spin waves [24] are propagating perturba- 
tions of the ordering in magnetic systems. These waves 
can have a variety of dispersion relations, depending on 
their particular type and propagation direction. The key 
points relevant to our present discussion are the follow- 
ing. First, several of these dispersion relations have the 
form u — ±(^1-1-^2^^), where Ai,A2 > 0; there is there- 
fore a lower frequency cutoff Ai. Second, a Doppler shift 
of spin waves can be induced by the application of an 
electric current [5S], an effect that was demonstrated a 
few years ago [26] . In the presence of an electric current 
the dispersion relation becomes w— Uk = ±{Ai + A2k'^), 
with a doppler-shifted frequency uj — Uk, where U is pro- 
portional to the current density [25j. One can therefore 
envisage creating horizons for spin waves using an electric 
current density that varies along the length of the sample 
due to a changing transverse area. Further details of this 
proposal will be given elsewhere. 



VIII. CONCLUSIONS 

The Hawking effect for massive Klein- Gordon modes 
does not exhibit a lower frequency cutoff given by the 
mass. This does not contradict the obvious fact that 
black holes cannot emit massive particles with uj > m 
to asymptotic inflnity. We have shown in detail that 
massive modes with u < m display a boomerang be- 
haviour: after being emitted from the horizon they even- 
tually turn around and are dragged back into the black 
hole. This bouncing, or blocking, behaviour is a familiar 
feature in analogue gravity in the context of waves with 
sub-luminal or super-luminal dispersion (or a combina- 
tion of both [T5]) at high wave numbers. In the Klein- 
Gordon case the interesting dispersive effects are at low 
wave numbers, but the same techniques of analysis reveal 
the behaviour. 

The lesson that the Hawking effect for Klein- Gordon 
modes is not suppressed compared to the massless case 
is an important one for analogue gravity. Hawking ra- 
diation of Klein-Gordon modes in an analogue system 
is just as amenable to experimental detection as radia- 
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tion of massless modes. Analogue-gravity systems with 
Klein-Gordon waves are possible using two-component 
BECs |21| . Langmuir waves in plasmas and barotropic 
waves in fluid mechanics. 
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